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We consider generic m x n-mode bipartitions of continuous variable systems, and study the associated bisym- 
metric multimode Gaussian states. They are defined as (m+n)-mode Gaussian states invariant under local mode 
permutations on the m-mode and n-mode subsystems. We prove that such states are equivalent, under local uni- 
tary transformations, to the tensor product of a two-mode state and of m + n — 2 uncorrected single-mode states. 
The entanglement between the m-mode and the n-mode blocks can then be completely concentrated on a single 
pair of modes by means of local unitary operations alone. This result allows to prove that the PPT (positivity 
of the partial transpose) condition is necessary and sufficient for the separability of (m + n)-mode bisymmet- 
ric Gaussian states. We determine exactly their negativity and identify a subset of bisymmetric states whose 
multimode entanglement of formation can be computed analytically. We consider explicit examples of pure and 
mixed bisymmetric states and study their entanglement scaling with the number of modes. 
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I. INTRODUCTION 

In quantum information and computation science, it is of 
particular relevance to provide theoretical methods to deter- 
mine the entanglement of systems susceptible to encompass 
many parties. Such an interest does not stem only from pure 
intellectual curiosity, but also from practical needs in the im- 
plementations of realistic information protocols. This is espe- 
cially true as soon as one needs to encode two-party informa- 
tion in a multipartite structure in order to minimize possible 
errors and decoherence effects (jl 0] • The study of the struc- 
ture of multipartite entanglement poses many formidable chal- 
lenges, concerning both its qualification and quantification, 
and so far little progress has been achieved for multi-qubit 
systems and in general for multi-party systems in finite dimen- 
sional Hilbert spaces. However, the situation looks somehow 
more promising in the arena of continuous variable systems, 
where some aspects of genuine multipartite entanglement can 
be at least qualitatively understood in the study of the entan- 
glement of multimode bipartitions. 

In the present work we study in detail the entanglement 
properties of multimode Gaussian states of continuous vari- 
able (CV) systems (for an introduction to CV quantum infor- 
mation see Ref. t3jj])- After the seminal analysis on the sep- 
arability of two-mode Gaussian states much progress 
has been accomplished on the separability conditions of mul- 
timode Gaussian states under various bipartitions 0000]. 
On the other hand, much less is known on the quantification of 
the entanglement of multimode, multipartite Gaussian states 
lHoll . In a previous work lUTll . we have presented a theoret- 
ical scheme to exactly determine the entanglement of, pure 
or mixed (n + l)-mode Gaussian states, under 1 x n-mode 
bipartitions, endowed with full or partial symmetries under 
mode exchange. More recently, a measure of genuine multi- 
partite CV entanglement has been proposed 11211 . that extends 
the approach introduced by Coffman, Kundu, and Wootters 
for multiqubit systems II 311 . and possesses a precise opera- 
tional meaning related to the optimal fidelity of teleportation 
in a continuous-variable teleportation network | l*4ll. 



In this paper we generalize the analysis introduced in 
Ref. Illll to bisymmetric (m + n)-mode Gaussian states of 
m x n-mode bipartitions. The main result of the present pa- 
per is that the bipartite entanglement of bisymmetric (m + n)- 
mode Gaussian states is unitarily localizable, i.e. that, through 
local unitary operations, it may be fully concentrated in a sin- 
gle pair of modes, each of them owned by one of the two par- 
ties (blocks). Here the notion of localizable entanglement is 
different from that introduced by Verstraete, Popp, and Cirac 
for spin systems [15]. There, it was defined as the maximal 
entanglement concentrable on two chosen spins through local 
measurements on all the other spins. Here, the local opera- 
tions that concentrate all the multimode entanglement on two 
modes are unitary and involve the two chosen modes as well, 
as parts of the respective blocks. 

The consequences of the unitary localizability are manifold. 
In particular, the PPT (positivity of the partial transpose) crite- 
rion is proved to be a necessary and sufficient condition for the 
separability of (m + n)-mode bisymmetric Gaussian states. 
Moreover, the block entanglement (i.e. the entanglement 
between blocks of modes) of bisymmetric (generally mixed) 
Gaussian states can be determined. The entanglement can be 
quantified by the logarithmic negativity in the general instance 
because the PPT criterion holds, but we will also show some 
explicit cases in which the entanglement of formation between 
m-mode and n-mode parties can be exactly computed. 

The plan of the paper is as follows. In Sec.|II]we introduce 
the notation and review some basic facts about Gaussian states 
and their entanglement properties. In Sec.|III]we show that a 
bisymmetric Gaussian state reduces to the tensor product of 
a correlated two-mode state and of uncorrelated single-mode 
states. In Sec. IIVI we exploit such a result to explicitly de- 
termine the entanglement of bisymmetric Gaussian states. In 
Sec.|V]the scaling of the block entanglement and the evalua- 
tion of the unitarily localizable entanglement involving differ- 
ent partitions of (generally mixed) symmetric states are stud- 
ied in detail. Finally, in Sec. IVII we present some conclusions 
and miscellaneous comments. 



2 



II. GAUSSIAN STATES OF BOSONIC SYSTEMS 

Let us consider a CV system, i.e. a system described by an 
infinite dimensional Hilbert space H = j Hi resulting 
from the tensor product of infinite dimensional Fock spaces 
Hi's. Let a, be the annihilation operator acting on Hi, and 
x% = (a; + a\) and pi = (cij — a\)/i be the related quadrature 
phase operators. The corresponding phase space variables will 
be denoted by Xi and pi. Let us group together the operators 
Xi and pi in a vector of operators X = (xi,pi, . . . , x n ,p n ). 
The canonical commutation relations (CCR) for the Xi's are 
encoded in the symplectic form ft 

[Xi,Xj] = 2iQij , 



with ft = u)® n , u = I _j q j . (1) 

A complete description of a CV quantum state g can be 
provided in terms of its symmetrically ordered characteris- 
tic function \- If we define the displacement operator D$ = 
exp(zJT T J7£), with £ S R 2 ™, then the characteristic function 
X associated to g is given by x(Q — Tr[g_D^]. The set of 
Gaussian states is, by definition, the set of states with Gaus- 
sian characteristic functions. Therefore a Gaussian state g is 
completely characterized by its first and second statistical mo- 
ments which form, respectively, the vector of first moments 

X = (jjti), (Xi), . . . , (X n ), {X n )J and the covariance ma- 
trix (CM) cr of elements 

o-y = i (XiXj + XjXi) - (Xi) (Xj) , (2) 

where, for any observable 6, (6) = Tr(gd), First statistical 
moments can be arbitrarily adjusted by local unitary opera- 
tions, which do not affect any property related to correlations 
or entropies. Therefore they will be unimportant to our aims 
and we will set them to in the following, without any loss of 
generality. Throughout the paper, cr will stand for the covari- 
ance matrix of the Gaussian state g. 

The positivity of g and the CCR entail the following relation 
on the CM cr of a quantum state g ("Robertson-Schrodinger" 
uncertainty relation) 

cr + ifl>0 7 (3) 

Inequality is the necessary and sufficient constraint cr has 
to fulfill to be a bona fide CM fl6ll . We mention that such a 
constraint implies cr > 0. 

The class of unitary transformations generated by sec- 
ond order polynomials in the field operators ('second-order' 
operations) is especially relevant in manipulating Gaussian 
states. For a n-mode systems, such operators may be mapped, 
through the so called 'metaplectic' representation, into the 
real symplectic group Sp(2 n ,ni) 11711 - made up by linear op- 
erations acting on a linear space (called 'phase space' in anal- 
ogy with classical Hamiltonian dynamics), which preserves 



the symplectic form ft under congruence: 

S € Sp(2n,W & S SIS = fl . 

Symplectic operations preserve the Gaussian character of the 
input state, acting linearly on first moments and by congru- 
ence on second moments: 

cr -» S T (TS . 

Ideal squeezers and beam splitters are examples of (respec- 
tively, 'active' and 'passive') symplectic transformations. 

A tensor product of Hilbert spaces (and of 'second-order' 
unitary operations) is mapped into a direct sum of phase 
spaces (and of symplectic transformations). Under a m X n 
mode partition, resulting from the direct sum of phase spaces 
Ti and T 2 with dimensions 2m and 2n respectively, we will 
refer to a transformation Si = Si © S2, with Si S Sp(2m,ffi) 
and S2 € Sp^nM) acting on Ti and T 2 , as to a "local sym- 
plectic operation". The corresponding unitary transformation 
is the "local unitary transformation" Ui = U\ <£> Ui. 

Let us recall that, due to a theorem by Williamson 1 the 
CM of a ?i-mode Gaussian state can always be written as fl^l 

cr = S T uS , (4) 

where S E Sp(2 n ,R) an d v is me CM 

v = diag(z/i,z/i, . . . ,v n ,v n ) , (5) 

corresponding to a tensor product of thermal states with diag- 
onal density matrix g® given by 




\k)i being the k-th number state of the Fock space Hi. The 
dual (Hilbert space) formulation of Eq. then reads: g = 
U* g® U, for some unitary U. The quantities i/,'s form the 
symplectic spectrum of the covariance matrix cr and can be 
computed as the eigenvalues of the matrix \iftcr\ \ 19]. Such 
eigenvalues are in fact invariant under the action of symplectic 
transformations on the matrix cr. 

The symplectic eigenvalues z/j encode essential informations 
on the Gaussian state g and provide powerful, simple ways 
to express its fundamental properties. For instance, provided 
that the CM cr satisfies cr > 0, then 

Ui > 1 

is equivalent to the uncertainty relation 0. We remark that 
the full saturation of the uncertainty principle can only be 
achieved by pure rt-mode Gaussian states, for which z/j = 
1 Vi = 1, . . . , n. Instead, mixed states such that Vi<k = 1 and 
Vi>k > 1, with 1 < fc < n, only partially saturate the uncer- 
tainty principle, with partial saturation becoming weaker with 
decreasing k. Such states are minimum uncertainty mixed 
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Gaussian states in the sense that the phase quadrature oper- 
ators of the first k modes satisfy the Heisenberg minimal un- 
certainty, while for the remaining n — k modes the state in- 
deed contains some additional thermal and/or Schrodinger- 
like correlations which are responsible for the global mixed- 
ness of the state. 

The symplectic eigenvalues are clearly invariant under sym- 
plectic operations. Yet, it is often advantageous to introduce 
other symplectic invariants, which can be easily handled in 
terms of second statistical moments. In the present work, deal- 
ing with a n-mode Gaussian state with CM er, we will make 
use of the obvious invariant Dct er (whose invariance is a con- 
sequence of the fact that Dct 5=1 V5 G Sp^n a)) an d of 
A„ = Ek =1 Det, 
of er: 



where the er^- are 2x2 submatrices 



O'nl 



(6) 



The invariance of A CT in the multimode case follows from its 
invariance in the case of two-mode states, proved in Ref. 1 20], 
and from the fact that any symplectic transformation can be 
decomposed as the product of two-mode transformations 12 ill . 
The symplectic eigenvalues of a two-mode Gaussian state 
are simply determined by the invariants introduced above: 



2(^) 2 = A CT T V&l -4 Deter 



(7) 



Also the purity fi = Tr g 2 of a multimode Gaussian state g, 
quantifying its degree of mixedness, is easily determined in 
terms of the symplectic invariants Det er, as 1 22] 



\i = 1/ v Dct er . 



(8) 



Regarding the entanglement of Gaussian states, we recall 
that the positivity of the partial transpose is a necessary and 
sufficient criterion for two-mode states to be separable (PPT 
criterion) |4]. The validity of such a criterion has been later 
extended to generic Gaussian states of 1 x n-mode systems 
0] and to (to + n)-mode Gaussian states with fully degen- 
erate symplectic spectrum I23t 12411 . For a bipartite system 
with Hilbert space H = Ha 8%, made up of two subsys- 
tems with Hilbert spaces Ha and Hb, the operation of partial 
transposition is defined as the transposition of the degrees of 
freedom associated to only one of the two subsystems, i.e. to 
the transposition of only one of the reduced Hilbert spaces, 
say Ha- Let us remark that the positivity of the partially 
transposed operator g does not depend on which subsystem 
is transposed nor on the basis chosen to perform the trans- 
position. Therefore the positivity of the partial transpose is 
invariant under local unitary transformations on the two sub- 
systems. In particular, for two-mode Gaussian states, the PPT 
criterion reduces to a simple inequality on the smallest sym- 
plectic eigenvalue v~ of the partially transposed CM er (par- 
tial transposition amounts to the mirror reflection of one of the 
four quadratures, see Ref. |4]). A two-mode Gaussian state is 
separable (i.e. not entangled) if and only if 



A proper quantification of the entanglement, easily com- 
putable for two-mode Gaussian states, is provided by the neg- 
ativity M , thoroughly discussed and extended in Ref. l25ll to 
CV systems (see also Refs. 12(1 127I0 . The negativity of a 
quantum state g is defined as 



YeWi - 1 



(10) 



where g is the partially transposed density matrix and ||6||i = 
Tr|6| stands for the trace norm of 6. The quantity Af(g) is 
equal to IX^-M' tne modulus of the sum of the negative 
eigenvalues of g, quantifying the extent to which g fails to 
be positive. Strictly related to M is the logarithmic negativity 
E_\f, defined as Em = ln||g||i, which constitutes an upper 
bound to the distillable entanglement of the quantum state q 
and is related to the entanglement cost under PPT preserving 
operations |28]. It can be easily shown 12911 that the loga- 
rithmic negativity of a two-mode Gaussian state is a simple 
function of the partially transposed symplectic eigenvalue i>~ 
alone: 



Em = max[0, — In v ] 



(11) 



quantifying the extent to which Inequality Q is violated. 

Let us recall that the bipartite entanglement of formation 
Ep l30ll of a quantum state q, shared by parties A and B, is 
defined as 

E f {q)= min yV#(|^)), (12) 

where the minimum is taken over all the pure states realiza- 
tions of g: 



and E(\4>i)) denotes the entropy of entanglement of the pure 
state | ?/>.;), defined as the von Neumann entropy of the reduced 
state obtained by tracing over the variables of one of the two 
subsystems: 

E{\^)) = -Tr A [Tr B |^)(^|ln(Tr B |^)(^|)]. 

As far as symmetric (i.e. with Det <j\\ = Det er 2 2, with ref- 
erence to the decomposition of Eq. ijS}) two-mode Gaussian 
states are concerned, the entanglement of formation Ep, can 
be computed 13111 . The quantity Ep turns out to be, again, a 
decreasing function of v~: 



Ef 



[0,h(D-)] , 



(13) 



with 



h(x) = £ + In f <-+->! in ( (1 - X? 



4.7.: 



4x 



4j: 



4.T 



> 1 



(9) 



Therefore the entanglement of formation provides, for two- 
mode symmetric Gaussian states, a quantification of entangle- 
ment fully equivalent to the one provided by the logarithmic 
negativity Ej\f. 
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III. STANDARD FORMS OF BISYMMETRIC 
MULTIMODE GAUSSIAN STATES 

We shall say that a multimode Gaussian state g is fully sym- 
metric if it is invariant under the exchange of any two modes. 
In the following, we will consider the fully symmetric m- 
mode and n-mode Gaussian states g a m and Qpn, with CMs 
(Tq,™ and crpn. Due to symmetry, we have that 



/a 



<T a m = 



\ 



c 



c 

V) 

(14) 

where a, e, (3 and £ are 2x2 real symmetric submatrices 
(the symmetry of e and £ stems again from the symmetry un- 
der the exchange of any two modes). All the properties related 
to correlations and entropic measures of multimode Gaussian 
states are invariant under local, single-mode symplectic op- 
erations. A first preliminary fact, analogous to the standard 
form reduction of two-mode states, will thus prove useful. 

Standard form of fully symmetric states. Let apn be the 

CM of a fully symmetric n-mode Gaussian state. The 2x2 
blocks f3 and C of er pn, defined by Eq. M4\ . can be brought 
by means of local, single-mode symplectic operations S € 
Spf^wi) i nt0 f orm ft = diag {b, b) and £ = diag [z\, Zi). 

Proof. The blocks (3, being CM's of reduced single mode 
Gaussian states, can be turned into their Williamson standard 
form by the same symplectic Si <G Sp(2,]R) acting on each 
mode. One is then left with the freedom of applying local, 
single-mode rotations that leave the blocks f3 invariant. The 
same rotation applied to each mode is sufficient to diagonalize 
C, since such a matrix is symmetric. □ 
The coefficients b, Z\, z 2 of the standard form are deter- 
mined by the local, single-mode invariant Det (5 = /ip 2 , 

and by the symplectic invariants Det <j pi = u^ 2 2 and A2 = 
A(crp2). Here up (/ipi) is the marginal purity of the single- 
mode (two-mode) reduced states, while A2 is the remain- 
ing invariant of the two-mode reduced states l32ll . This 
parametrization is provided, in the present instance, by the 
following equations 



1 

M/3 



21 



M/3 



(e_-e+), z 2 



VP 



(e_+e+), (15) 



with 



Mg2 



and 



\ 



M32 



This parametrization has a straightforward interpretation, be- 
cause fip and /j,p2 quantify the local mixednesses and A2 reg- 
ulates the entanglement of the two-mode blocks at fixed global 
and local purities l32ll . 



Let us next determine and analyse the symplectic spectrum 
(symplectic eigenvalues) of a pn. 

Symplectic degeneracy of fully symmetric states. The sym- 
plectic spectrum of cr pn is n — 1 times degenerate. The two 
symplectic eigenvalues of cr pn vZ and vt n read 



v = V '{b - zi)(b - z 2 ) 



J+ n = ^{b+{n-l)z 1 ){b+{n-l)z 2 ) , 
\)-times degenerate eigenvalue. 



(16) 



where Va is the (n - 

Proof. We recall that the symplectic eigenvalues of erpn are 
the absolute values of the eigenvalues of iflcrpn. Since the 
symplectic form f2 is block diagonal, with 2x2 blocks u> 
given by Eq. Q, the matrix iflcr is just the matrix er with iu> 
multiplying on the left any 2x2 block. Let us now consider 
the set of vectors {vi}, for i = 1, . . . , n — 1: 



(0,...,0,^, ^ ,0,...0) T 

modei modei+1 



(17) 



where, for covenience, we have introduced the two- 
dimensional vector v = (i b ~* 2 , 1) T . The Vi are n — 1 linear 

independent vectors. One has 

iSlavi = i(0, . . . ,0, H/3 - C)v) T , -(wG9 - C)v) T , 0, . . . , 0) T . 



mode i 



mode i+1 



(18) 



A straightforward computation gives 
iu>(f3-C)v = 




which recasts Eq. dl8l l into iflcrvi — v„ v.- L , thus proving that 
the symplectic eigenvalue of er is n — 1 times degener- 
ate. Note that, as one should expect, there exist also n — 1 
eigenvectors associated to the negative eigenvalue —vZ. To 

this end, it suffices to turn v into (— i _ 2 , 1) T . 

The remaining linearly independent eigenvector of iflcr pn 
is the vector 

(w T ,...,w T ) T , 



with w T — (iy/b+ (n — ijz\, y/b + (n — l)z 2 ) ■ 

It is immediate to verify that such a vector is associated to the 
eigenvalue vf n , completing the proof. □ 

The (n — l)-times degenerate eigenvalue v~Z is independent 
of n, while ut n can be simply expressed as a function of the 
single mode purity up and the symplectic spectrum of the two- 
mode block with eigenvalues v~p~ and v~p 2 '■ 



Hnf 



n(n - 2) (n- 1) / + 



n{v+ 2 y + {n-2)(vp) 2 ) . 

(20) 
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In turn, the two-mode symplectic eigenvalues are determined 
by the two-mode invariants by the relation 



2{vJ) 2 = Ap, T 



(21) 



The global purity Eq. (|8j of a fully symmetric multimode 
Gaussian state is 



(Det^»)- 1 / 2 = ((^)™- 1 4 



(22) 



and, through Eq. ( I20> . can be fully determined in terms of the 
one- and two-mode parameters alone. 

Obviously, analogous results hold for the m-mode CM 
er qtti of Eq. J14> . whose 2x2 submatrices can be brought 
to the form a, = diag (a, a) and e = diag (ei, e 2 ) and whose 
(to — 1) -times degenerate symplectic spectrum reads 



v a = (a-e 1 )(o-e 2 ) , 
y a'" = ( a + ( m _ l) e i)(a + (to — l)e 2 ) 



(23) 



Let us now generalize this analysis to the (to + n)-mode 
Gaussian states with CM er, which results from a correlated 
combination of the fully symmetric blocks er a m and er^: 



r 



(24) 



where T is a 2to x 2n real matrix formed by identical 2x2 
blocks 7. Clearly, T is responsible of the correlations exist- 
ing between the m-mode and the n-mode parties. Once again, 
the identity of the submatrices 7 is a consequence of the lo- 
cal invariance under mode exchange, internal to the m-mode 
and n-mode parties. States of the form of Eq. J24> will be 
henceforth referred to as bisymmetric . A significant insight 
into bisymmetric multimode Gaussian states can be gained by 
studying the symplectic spectrum of er and comparing it to the 
ones of a- a m and crpn. 

Symplectic degeneracy of bisymmetric states. The symplec- 
tic spectrum of the CM er Eq. \24\ of a bisymmetric (m + n)- 
mode Gaussian state includes two degenerate eigenvalues, 
with multiplicities to — 1 and n — 1. Such eigenvalues co- 
incide, respectively, with the degenerate eigenvalue v~ of the 
reduced CM a a m and the degenerate eigenvalue v7 of the 
reduced CM (Tpm. 

Proof. One can proceed constructively, in analogy with the 
proof of the previous proposition. Let us consider the standard 
forms of the blocks cr a m and crpn, while keeping the 2x2 
submatrices 7 in arbitrary, generally nonsymmetric, form. Let 
us next focus on the block crpn and define the vectors Vi by 



(o,...,o,iO 



(25) 



They are the vectors obtained from the vectors Vi's of Eq. illi 
by appending to them 2to null entries on the left. Because of 
the identity of the blocks 7, their contributions to the secu- 
lar equation cancel out and it is straightforward to verify that 
the vectors Vi's are n — 1 eigenvectors of ifier with eigen- 
value vZ . The same argument holds considering the subma- 
trix cr Qm , thus completing the proof. □ 



Equipped with these results, we are now in a position to de- 
termine the bipartite entanglement of bisymmetric multimode 
Gaussian states and prove that it can always be unitarily local- 
ized or concentrated. 

Unitary localization of the entanglement of bisymmetric 
states. The bisymmetric (m + n)-mode Gaussian state with 
CM er Eq. \24\ can be brought, by means of a local unitary 
operation, with respect to the to x n-mode bipartition with 
reduced CMs cr a m and crpn, to a tensor product of single- 
mode uncorrelated states and of a two-mode Gaussian state. 

Proof. Let us focus on the n-mode block crpn . The vectors Vi 
of Eq. d25i . with the first 2m entries equal to 0, are, by con- 
struction, simultaneous eigenvectors of iftcpn and iJler, with 
the same (degenerate) eigenvalue. This fact suggests that the 
phase-space modes corresponding to such eigenvectors are the 
same for er and for crpn. Then, bringing by means of a local 
symplectic operation the CM trpn in Williamson form, any 
(2n — 2) x (2n — 2) submatrix of er will be diagonalized be- 
cause the normal modes are common to the global and local 
CMs. In other words, no correlations between the m-mode 
party with reduced CM cr a m and such modes will be left: all 
the correlations between the m-mode and n-mode parties will 
be concentrated in the two conjugate quadratures of a single 
mode of the n-mode block. Going through the same argu- 
ment for the m— mode block with CM er Qm would prove the 
proposition and show that the whole entanglement between 
the two multimode blocks can always be concentrated in only 
two modes, one for each of the two multimode parties. 

To prove this property we proceed first by investigating the 
relationship between the transformations which diagonalize 
iftcr and the symplectic operations that bring er in Williamson 
normal form u 1 33 ] . The problem one is immediately faced 
with is that these transformations are not unique because the 
normal form associated to er is invariant under local rotations 
(this local freedom is always present in the selection of normal 
modes) and, due to degeneracy, also under global symplectic 
rotations of the modes associated to the degenerate eigenvalue 

. Thus there is an ambiguity in selecting the eigenvectors 
of ifler and therefore in determining the transformation that 
diagonalizes it. Moreover, if {wi} is a set of 2 (m+n) column- 
vectors normalized eigenvectors of zfier, then any matrix T of 
the form 



(26) 



diagonalizes iflcr: T -1 (if2er)T = D (with the £j's arbitrary 
complex coefficients). However, we can proceed by observ- 
ing that the 2x2 matrix ia> is diagonalized by the unitary 
transformation U, with 



u= — ( i ~ l 

M 1 1 



so that UHaujU — diag (a, —a) (where a is any complex 
number). We can then define the matrix U — fj®( m + n ) j 
which is local in the sense that it is block diagonal and acts 
on each mode separately, such that for any normal form v 



U-HftvU = D. 



(27) 
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where D = T^ 1 iflcrT is a diagonal matrix with entries 
{T^i} ( m terms of the symplectic eigenvalues). Let us next 
denote by S one of the symplectic transformations that bring 
<t in normal form: S T crS = v. It is then easy to see that 



D = T~ l {in(T)T = U^iiflup 

= U- 1 {iQ.S T aS)U = U^S^iinajSU 



and therefore 



S = TU- 



(28) 



(29) 



where in Eq. (12 8 1 we have exploited the fundamental property 
of symplectic transformations: S^^ttS^ 1 = CI. Eq. d29l 
shows that there must exist some symplectic transformation 
that diagonalizes iClcr and satisfies the further condition given 
by Eq. d29i . In fact, it is obvious that not every T diagonal- 
izing itlcr is a symplectic transformation when multiplied on 
the right by f/T. Viceversa, if this last condition holds, the 
symplectic operation that brings er in normal form is given by 
Eq. (I29> . The modes that diagonalize the quadratic form er in 
phase space can be reconstructed in terms of S: since they are 
linear combinations of the original modes and S T crS is diag- 
onal, they can be expressed by real column vectors identified 
by the columns of S. 

We can now go back to our original problem: leaving aside 
the involved task of exactly determining which choice of the 
eigenvectors of ifier leads to a symplectic transformation of 
the form Eq. (1291 . we are anyway assured that in the sub- 
space associated to the eigenvalues ^vZ such eigenvectors 
must be linear combinations of the Wj's defined in Eq. (125 \ 
and their counterparts associated to the eigenvalue — vZ (with 
their first 2m entries, related to the rn-mode party, set equal to 
0). Therefore the transformation T reads, in general, 



T 



1,771 



T m , 1 

T m +l.l 



m+n,l 



T 



o 
o 

m+l,m+l 



o 



T 



m-\-n,m ^m+n,m+l 



-^m+l,m-f n 

(30) 

where stands for 2x2 null matrices and Tj j are 2x2 
blocks, whose exact form is unessential to our aims. Exploit- 
ing Eq. (1291 . for the last 2(n — 1) columns of S we obtain, in 
terms of 2 x 2 matrices, 



t \T 



( 0___0 ,U*Tl t ,...,U*Tl t ) 

first m modes 



(31) 



Due to the presence of the first m null entries, the n — 1 modes 
determined by Eq. J3 1 1 are normal modes of both the global 
CM <j and the local CM trpn. An analogous proof, going 
along the same lines of reasoning, holds for the reduced CM 
cr a m : it can be reduced to a local normal form that shares m — 
1 normal modes with the global CM er. These results imply 
that the form in which all the correlations between the two 
parties are shared only by a single mode of the n-mode party 



and by a single mode of the m-mode party can be obtained by 
means of local symplectic (unitary) operations, namely by the 
symplectic operations bringing the block crpn and the block 
cr am in Williamson form. 

For ease of the reader and sake of pictorial clarity, we can 
supplement the proof by explicitly writing down the different 
forms of the CM <r at each step; such matrix representations 
allow an immediate visualization of the process of unitary 
concentration of the entanglement between a single pair of 
modes, one for each multimode party. The CM a of a bisym- 
metric (m + n)-mode Gaussian state reads (see Eq. d24l ) 
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According to what we have just shown, reducing to normal 
form the block crpn brings the global CM a in the form CM 
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where the 2x2 blocks i/T„ = vf n I2 and u~Z = vZ I2 are 
the Williamson normal blocks associated to the two symplec- 
tic eigenvalues of apn. The identity of the submatrices 7' is 
due to the invariance under permutation of the first m modes, 
which are left unaffected. The subsequent symplectic diago- 
nalization of cr a m puts the global CM a in the following form 
(notice that the first, m + 1-mode reduced CM is again a ma- 
trix of the same form of er, with n = 1): 
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with i/Jm = i/^n 1 2 and i/~ = i/ l 2 . Eq. (l33l shows explic- 
itly that the state with CM cr", obtained from the original state 
with CM cr by exploiting local unitary operations, is the ten- 
sor product of to + n — 2 uncorrelated single-mode states and 
of a correlated two-mode Gaussian state. The proof is there- 
fore complete, and shows that the amount of entanglement 
(quantum correlations) present in any bisymmetric multimode 
Gaussian state can be localized (concentrated) in a two-mode 
Gaussian state (i.e. shared only by a single pair of modes), via 
local unitary operations. These results and their consequences 
will be discussed in detail in the following sections. □ 



IV. BLOCK ENTANGLEMENT OF MULTIMODE 
GAUSSIAN STATES 

In the previous section, the study of the multimode CM cr of 
Eq. (I32> has been reduced to a two-mode problem by means of 
local unitary operations. This finding allows for an exhaustive 
analysis of the bipartite entanglement between the m- and n- 
mode blocks of a multimode Gaussian state, resorting to the 
powerful results available for two-mode Gaussian states. For 
any multimode Gaussian state with CM cr, let us define the 
associated equivalent two-mode Gaussian state g eq , with CM 
(T eq given by 

= ( i ) ' (34) 

where the 2x2 blocks have been implicitly defined in the CM 
( 1331 . As already mentioned, the entanglement of the bisym- 
metric state with CM cr, originally shared among all the m+n 
modes, can be completely concentrated by local unitary (sym- 
plectic) operations on a single pair of modes in the state with 
CM cr eq . Such an entanglement is, in this sense, localizable. 
Obviously, this kind of localization of entanglement by local 
unitaries is conceptually very different from the localization 
of entanglement by local measurements first discussed by Ver- 
straete, Popp, and Cirac for qubit systems Hal . We now move 
on to describe some consequences of this result. 

A first qualificative remark is in order. It is known that the 
PPT criterion is necessary and sufficient for the separability of 
Gaussian states of 1 x 1-mode and 1 x n-mode bipartitions. In 
view of the invariance of such a criterion under local unitary 
transformations, which can be appreciated by the definition 
of partial transpose at the Hilbert space level, and considering 
the results proved in the previous section, it is immediate to 
verify that the following property holds: 

PPT criterion for bisymmetric multimode Gaussian states. 

For generic m x n-mode bipartitions, the positivity of the par- 
tial transpose (PPT) is a necessary and sufficient condition for 
the separability of bisymmetric m + n-mode Gaussian states. 

This statement is a first important generalization to to x n 
bipartitions of the result proved by Werner and Wolf for the 
case of 1 x n bipartitions [6]. In particular, it implies that 
no bisymmetric bound entangled Gaussian states may exist 



and all the m x n block entanglement of such states 
is distillable. Moreover, it justifies the use of the negativity 
and the logarithmic negativity as measures of entanglement 
for these multimode Gaussian states. 

As for the quantification of the entanglement, exploiting 
some recent results on two-mode Gaussian states 11291 13211 we 
can select the relevant quantities that, by determining the cor- 
relation properties of the two-mode Gaussian state with CM 
cr eq , also determine the entanglement and correlations of the 
multimode Gaussian state with CM cr. These quantities are, 
clearly, the equivalent marginal purities [i aeq and np eq , the 
global purity u eq and the equivalent two-mode invariant A eq . 
Let us remind that, by exploiting Eqs. (I16> . J23i and (I15> . the 
symplectic spectra of the CMs cr a m and cr^n may be recov- 
ered by means of the local two-mode invariants up, fi a , [ia2, 
fi a 2, Afj2 and A a 2. The quantities /.i aeq and p,p eq are easily 
determined in terms of local invariants alone: 

Mae 9 = Up eq = l/f£ n ■ (35) 

On the other hand, the determination of fi eq and A eq require 
the additional knowledge of two global symplectic invariants 
of the CM cr; this should be expected, because they are sus- 
ceptible of quantifying the correlations between the two par- 
ties. The natural choices for the global invariants are the 
global purity fx = 1/ VDct cr and the invariant A, given by 

A = wDet a + m(m — l)Det e + nDet j3 
+n(n — l)Det £ + 2mnDetf ■ 

One has 

^ = o^r-^r-v , (36) 

A eq = A - (to - l)(v-) 2 - (n - 1)(^) 2 . (37) 

The entanglement, quantified by the logarithmic negativity, 
and the mutual information between the TO-mode and the n- 
mode subsystems can thus be easily determined, as it is the 
case for two-mode states. In particular, the smallest sym- 
plectic eigenvalue v eq of the matrix er eg , derived from or eq 
by partial transposition, fully quantifies the entanglement be- 
tween the TO-mode and n-mode partitions. Recalling the re- 
sults known for two-mode states l29l l32ll . the quantity v eq 
reads 

2D 2 - A - A 2 - — 

V t^eq 

2 2 

with A eq = — h — A eq . 

l^aeq l^fjeq 

The logarithmic negativity measuring the bipartite 

entanglement between the m-mode and n-mode subsystems 
is then 

Efj- = max [— In v eq , 0] . (38) 

In the case v^m = ut„ , corresponding to the condition 

(o+(m-l)e 1 )(a+(m-l)e 2 ) = (b+(n-l)z 1 )(b+(n~l)z 2 ) , 

(39) 
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the equivalent two-mode state is symmetric and we can de- 
termine also the entanglement of formation, using Eq. dl 3I >. 
Let us note that the possibility of exactly determining the en- 
tanglement of formation of a multimode Gaussian state of a 
m x n-mode repartition is a rather remarkable consequence, 
even under the symmetry constraints obeyed by the CM <x. 
Another relevant fact to point out is that, since both the loga- 
rithmic negativity and the entanglement of formation are de- 
creasing functions of the quantity v eq , the two measures in- 
duce the same entanglement hierarchy on such a subset of 
equivalently symmetric states {i.e. states whose equivalent 
two-mode CM cr eq is symmetric). 



From Eq. d36i it follows that, if the (m + n)-mode bisym- 



metric state is pure (p 



1), then the equiva- 



lent two-mode state is pure as well (p eq = 1) and, up to lo- 
cal symplectic operations, it is a two-mode squeezed vacuum. 
Therefore any pure bisymmetric multimode Gaussian state is 
equivalent, under local unitary (symplectic) operations, to a 
tensor product of a pure two-mode squeezed vacuum and of 
m + n — 2 uncorrelated vacua. 



More generally, if both the reduced m-mode and ?i-mode 
CMs (T a m and crpm of a bisymmetric, mixed multimode 
Gaussian state cr of the form Eq. 1241 correspond to Gaussian 
mixed states of partial minimum uncertainty, i.e. if v~ v 



1, then Eq. (1361 implies p eq = p. Therefore, the 
equivalent two-mode state has the same entanglement and the 
same degree of mixedness of the original multimode state. In 
all other cases of bisymmetric multimode states one has that 
p eq > p and the process of localization produces a two- 
mode state with higher purity than the original multimode 
state. In this specific sense, we see that the process of lo- 
calization implies a process of purification as well. We can 
understand this key point observing that the entanglement is 
localized by performing local unitary transformation which 
are reversible by definition. Then, in principle, by only using 
passive and active linear optics elements such as beam split- 
ters, phase shifters and squeezers |9], one can implement a 
reversible machine that, from mixed, bisymmetric multimode 
states with strong quantum correlations between all the modes 
(and consequently between the m-mode and the n-mode par- 
tial blocks) but weak couplewise entanglement, is able to ex- 
tract a highly pure, highly entangled two-mode state (with no 
entanglement lost, all the m x n entanglement can be local- 
ized). If needed, the same machine would be able, starting 
from a two-mode squeezed state and a collection of uncorre- 
lated thermal or squeezed states, to distribute the two-mode 
entanglement between all modes, converting the two-mode 
into multimode, multipartite quantum correlations, again with 
no loss of entanglement. The bipartite or multipartite entan- 
glement can then be used on demand, the first for instance in a 
CV quantum teleportation protocol, the latter to secure quan- 
tum key distribution or to perform multimode entanglement 
swapping. 



V. QUANTITATIVE LOCALIZATION OF THE BLOCK 
ENTANGLEMENT 

In this section we will explicitly compute the block en- 
tanglement (i.e. the entanglement between different blocks of 
modes) for some instances of multimode Gaussian states. We 
will study its scaling behavior as a function of the number 
of modes and explore in deeper detail the localizability of 
the multimode entanglement. We focus our attention on fully 
symmetric 2n-mode Gaussian states described by a 2n x 2n 
CM <j pin given by Eq. J14t . These states are trivially bisym- 
metric under any bipartition of the modes, so that their block 
entanglement is always localizable by means of local sym- 
plectic operations. Let us recall that concerning the covari- 
ances in normal forms of fully symmetric states (see Sec. Will , 
pure states are characterized by 

Zi = [l + 6 2 (27i-2)-(2n-l)-(-l) J 



x v /(6 2 -l)((26n) 2 -(2n-2) 2 ) /[26(2n-l)] , 

(40) 

and belong to the class of CV GHZ-type states discussed in 
Refs. I9l ll ill . These multipartite entangled states are gener- 
ated as the outputs of the application of a sequence of 2n — 1 
beam splitters to 2n single-mode squeezed inputs |9]. In the 
limit of infinite squeezing, these states reduce to the simulta- 
neous eigenstates of the relative positions and the total mo- 
mentum, which define the proper GHZ states of CV systems 
J9J]. The CM er^ 2 „ of this class of pure states, for a given num- 
ber of modes, depends only on the parameter b = 1/ pp > 1, 
which is an increasing function of the single-mode squeezing. 
Correlations between the modes are induced according to the 
above expression for the covariances Zj. Exploiting our pre- 
vious analysis, we can compute the entanglement between a 
block of k modes and the remaining 2n — k modes, both for 
pure states (in this case the block entanglement is simply the 
Von Neumann entropy of each of the reduced blocks) and, re- 
markably, also for mixed states. 

We can in fact consider a generic 2n-mode fully symmetric 
mixed state with CM <x^, obtained from a pure fully sym- 
metric (2n + g)-mode state by tracing out q modes. For any 
q, for any dimension k of the block (k < ri), and for any non 
zero squeezing (i.e. for b > 1) one has that V). < 1, mean- 
ing that the state exhibits genuine multipartite entanglement, 
as first remarked in Ref. |9] for pure states: each fc-mode 
party is entangled with the remaining (2n — fc)-mode block. 
Furthermore, the genuine multipartite nature of the entangle- 

0k I o2n-fc 

ment can be precisely quantified by observing that Eft- 
is an increasing function of the integer k < n, as shown in 
Fig.^ Moreover, we note that the multimode entanglement of 
mixed states remains finite also in the limit of infinite squeez- 
ing, while the multimode entanglement of pure states diverges 
with respect to any bipartition, as shown in Fig.^ 

In fully symmetric Gaussian states, the block entanglement 
is localizable with respect to any k x (2n — k) bipartition. 
Since in this instance all the entanglement can be concentrated 
on a single pair of modes, after the partition has been decided, 
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FIG. 1: (color online). Hierarchy of block entanglements of fully 
symmetric 2n-mode Gaussian states of k x (2n — k) bipartitions 
(n — 10) as a function of the single-mode squeezing b. The block 
entanglements are depicted both for pure states (solid lines) and for 
mixed states obtained from fully symmetric (2n + 4) -mode pure 
Gaussian states by tracing out 4 modes (dashed lines). All the quan- 
tities plotted are dimensionless. 



no strategy could grant a better yield than the local symplectic 
operations bringing the reduced CMs in Williamson form (be- 
cause of the monotonicity of the entanglement under general 
LOCC). However, the amount of block entanglement, which 
is the amount of concentrated two-mode entanglement after 
unitary localization hase taken place, actually depends on the 
choice of a particular k x (2n — k) bipartition, giving rise to 
a hierarchy of localizable entanglements. 

Let us suppose that a given Gaussian multimode state (say, 
for simplicity, a fully symmetric state) is available and its en- 
tanglement is meant to serve as a resource for a given proto- 
col. Let us further suppose that the protocol is optimally im- 
plemented if the entanglement is concentrated between only 
two modes of the global systems, as it is the case, e.g., in a 
CV teleportation protocol between two single-mode parties. 
Which choice of the bipartition between the modes allows for 
the best entanglement concentration by a succession of local 
unitary operations? In this framework, for an even number 
of modes, the worst localization strategy consists in assign- 
ing k = 1 mode at one party and 2n — 1 modes to the other. 
Conversely, the best option for localization is an equal k = n 
splitting of the In modes between the two parties. The log- 
arithmic negativity , concentrated into two modes by 
local operations, represents the optimal localizable entangle- 
ment (OLE) of the state trasn, where "optimal" refers to the 
choice of the bipartition. Clearly, the OLE of a state with 

In + 1 modes is given by EKr 1 . These results may be 
applied to arbitrary, pure or mixed, fully symmetric Gaussian 
states. 

We now turn to the study of the scaling behavior with n of 
the OLE of 2n-mode states, to understand how the number 
of local cooperating parties can improve the maximal entan- 
glement that can be shared between two parties. For generic 
(mixed) fully symmetric 2n-mode states ofnxn bipartitions, 
the OLE can be quantified also by the entanglement of forma- 
tion Ep, as the equivalent two-mode state is symmetric. It is 
then useful to compare, as a function of n, the lxl entan- 



FIG. 2: (color online). Scaling, with half the number of modes, of 
the entanglement of formation in two families of fully symmetric 
2n-mode Gaussian states. Diamonds denote pure states, while mixed 
states (denoted by stars) are obtained from (2n + 4)-mode pure states 
by tracing out 4 modes. For each class of states, two sets of points 
are plotted, one referring tonxn entanglement (filled symbols), and 
the other to 1 x 1 entanglement (empty symbols). Notice how the 
n x n entanglement, equal to the optimal localizable entanglement 
(OLE) and estimator of genuine multipartite quantum correlations 
among all the 2n modes, increases at the detriment of the bipartite 
lxl entanglement between any pair of modes. The single-mode 
squeezing parameter is fixed at b = 1.5. All the quantities plotted 
are dimensionless. 



glement of formation between a pair of modes (all pairs are 
equivalent due to the global symmetry of the state) before the 
localization, and the n x n entanglement of formation, which 
is equal to the optimal entanglement concentrated in a specific 
pair of modes after performing the local unitary operations. 
The results of this study are shown in Fig. [2] The two quanti- 
ties are plotted at fixed squeezing b as a function of n both for 
a pure 2n-mode state with CM <Xo 2ra and a mixed 2n-mode 



state with CM tr^t ■ As the number of modes increases, any 
pair of modes becomes steadily less entangled, but the total 
multimode entanglement of the state grows and, as a conse- 
quence, the OLE increases with n. In the limit n — > oo, the 
n x n entanglement diverges while the lxl one vanishes. 
This holds both for pure and mixed states, although the global 
degree of mixedness produces the typical behavior that tends 
to reduce the total entanglement of the state. 



VI. CONCLUDING REMARKS 

We have shown that bisymmetric multimode Gaussian 
states (pure or mixed) can be reduced, by local symplectic 
operations, to the tensor product of a correlated two-mode 
Gaussian state and of uncorrelated thermal states (the latter 
being obviously irrelevant as far as the correlation proper- 
ties of the multimode Gaussian state are concerned). As a 
consequence, all the entanglement of bisymmetric multimode 
Gaussian states of arbitrary m x n bipartitions is unitarily lo- 
calizable in a single (arbitrary) pair of modes shared by the 
two parties. Such a useful reduction to two-mode Gaussian 
states is somehow similar to the one holding for states with 
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fully degenerate symplectic spectra l23l l24ll . encompassing 
the relevant instance of pure states, for which all the symplec- 
tic eigenvalues are equal to 1, The present result allows to 
extend the PPT criterion as a necessary and sufficient condi- 
tion for separability for all bisymmetric multimode Gaussian 
states of arbitrary m X n bipartitions, and to quantify their 
entanglement. 

Notice that, in the general bisymmetric instance addressed 
in this work, the possibility of performing a two-mode reduc- 
tion is crucially partition-dependent. However, as we have 
explicitly shown, in the case of fully symmetric states all the 
possible bipartitions can be analysed and compared, yielding 
remarkable insight into the structure of the multimode block 
entanglement of Gaussian states. This leads finally to the de- 
termination of the maximum, or optimal localizable entangle- 
ment that can be concentrated on a single pair of modes. 

It is important to notice that the multipartite entanglement 



in the considered class of multimode Gaussian states can be 
produced and detected |S| 13511 . and also, by virtue of the 
present analysis, reversibly localized by all-optical means. 
Moreover, the multipartite entanglement allows for a reliable 
{i.e. with fidelity T > T c , where T c = 1/2 is the classi- 
cal threshold) quantum teleportation between any two parties 
with the assistance of the remaining others 

El- 

This quantum 

teleportation network has been recently demonstrated experi- 
mentally with the use of fully symmetric three-mode Gaussian 
states 13611 - 
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